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Abstract 

In this work we focus on various phenomenological aspects of the lightest even tiers, 
(2,0) and (0,2), in models based on a Real Projective Plane in 6 dimensions. We 
discuss the spectrum of the levels due to loop corrections, and the limit when the two 
radii are equal, in which case the two levels mix with each other and a new basis is 
defined. We also discuss the dependence of the spectrum on the ratio of the two radii. 
These results are essential to understand the phenomenology of the model at colliders 
(LHC) and to predict the relic abundance of Dark Matter. Finally, we estimate the 
bounds on the radius from resonant decays of the even tiers at the LHC, showing that 
they can be in the 600 GeV range after the complete analysis of the 2011 data. 



1 Introduction 



We live in exciting times for particle physics! Since its first physics run, the Large Hadron 
Collider (LHC) experiments have been collecting a considerable amount of high quality data. 
The results so far have been disappointing for Beyond the Standard Model (BSM) physics, 
because the first searches have accurately confirmed the expectations for a Standard Model 
(SM) only hypothesis, and no trace of new physics have appeared. Many models, like simple 
or constrained supersymmetric models, Randall-Sundrum gravitons, large extra dimensions 
(d la ADD), TeV-scale Black Holes, excited fermions and leptoquarks, fourth generation 
of fermions, have been disfavoured or severely constrained. However, this was the easy 
catch, selected for an early discovery, and many simple and motivated models are still poorly 
constrained by the present data. The ATLAS and CMS collaborations have also announced 
the discovery of a new resonance at a mass of 125 -v- 126 GeV, that has similar properties as 
the SM Higgs boson, thus confirming the hints present in the 2011 dataset. Some interesting 
discrepancies are present in the measured rates of various channels, mainly photons and tau 
ones, however more statistics will be necessary in order to derive any clear indication on the 
nature of the resonance. However, even if the discovery were finally associated with the SM 
Higgs and the absence of new physics signals were to persist, the interest in extensions of 
the SM would not fade. In fact, it will take much more effort to determine if the resonance 
discovered by the experiments is the SM Higgs or not. Furthermore, the discovery of new 
particles in the messy hadronic environment at the LHC experiments can be more challenging 
than we hoped for: new physics may well require more statistics or more sophisticated search 
strategies. 

A compelling evidence of new physics has been around since 80 years: it is now confirmed 
in a wide range of observations, from astrophysical observations of the rotation curves in 
galaxies, through weak gravitational lensing and the Cosmic Microwave Background fluc- 
tuations, up to simulations of the galaxy formations in the early Universe, that we need a 
large amount of Dark Matter in the Universe. While a new paradigm for our understanding 
of the dynamics of the Universe may be required, the puzzle can be easily solved by the 
presence of a stable massive new particle, unaccounted for in the SM, which, besides gravi- 
tation, interacts weakly with standard matter. Motivated by this observation, many models 
of new physics have been proposed which have a Dark Matter candidate: a rich zoology of 
such reclusive beasts can be found in the literature. A common ground is the presence of 
a discrete symmetry that stabilises the Dark Matter particle and prevents or strongly sup- 
presses decays to SM states: such symmetry is usually added ad hoc, and in the best cases is 
accidental or it carries along other benefits for the model. For instance, in supersymmetry, 
R-parity can also prevent baryon number violating terms in the supersymmetric Lagrangian, 
therefore forbidding proton decay. Extra dimensions have a natural set of symmetries that 
may play the role of the Dark Matter parity: namely symmetries of the compact space, 
or Kaluza-Klein (KK) parities jl]. However, it has been shown in j2] that most scenarios 
studied in the literature require additional conditions on the 4-dimensional counter-terms lo- 
calised on the boundaries or singular points of the space, therefore the parity is not a direct 
consequence of the extra dimensional set-up and should be generically considered ad-hoc. 
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Furthermore, the KK-parity is easily broken once interesting models are constructed: for 
instance, Gauge-Higgs unification models require the localisation of fermion zero modes on 
4D fixed points [3], thus potentially breaking the symmetry. 

For instance, in 5 dimensional models compactified on an interval 5' 1 /Z 2 , the KK parity 
is given by a translation by the length of the interval (frequently defined as an inversion 
with respect to the centre of the interval, which is an equivalent definition but which do not 
act properly on the fields). The introduction of a bulk mass term for fermions, which is odd 
under the Z 2 that defines the orbifold and would localise the chiral zero mode on either end 
of the interval depending on the sign of the mass term, has been shown to be incompatible 
with the inversion symmetry 14]. The main issue seems to be related with the presence 
of fixed points or boundaries, which are however a common feature to all compactifications 
with chiral zero modes. This drawback can be turned into a powerful tool by the observation 
that the requirement of an exact KK-parity is very restrictive and can be used as a powerful 
selection criterion to single out promising compact spaces. This criterion has been used in [2] 
to show that there exists only one 6D orbifold in 5 or 6 flat dimensions that possesses such 
a symmetry: the real projective plane (RP 2 ). 

Some aspects of the phenomenology of the RP 2 have been discussed in (2j and [5]. In 
this paper we will focus on the even tiers (2,0) and (0,2): they are important because they 
lead to events without missing energy, and therefore they can give signatures very different 
from supersymmetric models. Furthermore, they may play a crucial role in this model where 
naturally small mass splittings are generated by loop corrections, thus the phenomenology of 
the lightest odd tier resembles a compressed spectrum supersymmetric scenario [6j|7] . Due to 
the large interest given to supersymmetric models (both ATLAS and CMS have a dedicated 
team on supersymmetric, while all other models are binned in the Exotica group), many 
of these signatures have not been thoroughly explored. In this paper, we will discuss for 
the first time the effect of the two different radii allowed by the compactification and the 
structure of the symmetries of the space. Then we will introduce the results for the one loop 
corrections to the masses and discuss in details the spectrum as a function of the two radii. 
This discussion is crucial in the understanding of the LHC phenomenology of the model. 

The paper is organised as follows: after introducing the model in Section [2j we present a 
detailed discussion of the mass spectrum of the tiers (2,0) and (0,2) in Section[3j In Section|i] 
we discuss how we implemented the model in Monte Carlo generators, and in Section [5] we 
show a series of estimates of the bounds from present LHC searches before concluding in 
Section [6j 

2 The real projective plane 

We will consider a quantum field theory defined on a d dimensional flat manifold chosen 
to be the direct product of the standard four-dimensional Minkowski space-time and a 
d — 4-dimensional orbifolcfj] Our aim is to consider orbifolds without fixed points. Fixed 
points are not dangerous by themselves, however, as in general they break the <i-dimensional 

*An orbifold is defined as a quotient space of a manifold modulo a discrete symmetry group 
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Lorentz invariance down to the 4-dimensional one, the predictivity of the theory defined 
on an orbifold with fixed points will be limited. This limitation arises from the fact that 
the divergences appearing when calculating loop corrections require counter-terms localised 
on these points. Furthermore, any symmetry associated with the geometry of the space, 
including the eventual KK-parity, will be broken by generic localised terms unless one imposes 
some ad hoc symmetry conditions relating the singularities. As the extra dimensional theories 
are appealing mainly because they provide a stable dark matter candidate we want the 
symmetry preserving the stability of the dark matter particle to be an inherent property of 
the space. 

The simplest case would be to study the 1-dimensional orbifolds, that is a circle or an 
interval. The circle is the simplest orbifold that can be defined as S 1 = M./Z, that is by 
identifying the points y ~ y + 2ttR on an infinite line 1R (R is the radius of the circle and y 
the extra space co-ordinate). It is known that this orbifold has no fixed points, but no chiral 
fermion can be defined on it. Therefore this possibility is ruled out. Chiral fermions can be 
defined on the interval S 1 /^, where the Z2 symmetry identifies points y ~ — y = 2n — y: 
this orbifold, however, has two fixed points at the boundaries of the interval, y = and 
y = it R. From a geometrical point of view, the two fixed points are in-equivalent: they are 
defects of the space (or branes) where localised 4-dimensional interactions can be added. 
Such interactions, which are required as counter-terms for logarithmic divergences at one 
loop, do violate all the symmetries of the extra dimension, therefore no KK parity will 
in general survive. The bulk of the interval is invariant under a mirror symmetry around 
the centre of the interval, y — > ttR — y, which however interchanges the two fixed points. 
Therefore, imposing by hand that the localised interactions on the two fixed points are equal, 
a KK parity can be obtained: the mirror symmetry is not a good symmetry because the 
two chiralities of massive fermions would have opposite parity, however it is equivalent to 
a translation y — >■ y + tcR combined with the Z2 symmetry that defines the orbifold. This 
scenario has been studied in the literature |8j , and it was the first example of Dark Matter 
in extra dimensions: however, the symmetry is imposed ad hoc and it is easily broken, for 
instance by bulk mass terms for the fermions which control the localisation of the massless 
zero modes. 

The other face of the coin is that the requirement of spaces having an exact KK parity is 
a very effective selection rule on the number of viable orbifolds. As it was pointed out in [2] 
there is a unique 2-dimensional orbifold, among the 17 in-equivalent orbifolds that can be 
defined on a 2-dimensional euclidean plane, that has no fixed points and such that massless 
chiral fermionic fields can be defined on it: the real projective plane. The real projective 
plane is a compact, non-orientable orbifold of Euler characteristic 1 without boundaries. 
Among all possible descriptions of the real projective plane that are topologically equivalent 
spaces, we choose the simplest case with flat bulk metric: the spherical projective plane has 
been considered in [9], and it leads to a completely different phenomenology. 
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Figure 1: The fundamental domain of the flat real projective plane of side lengths of 7ri?4 and 7ri?5 (within 
black thick dashed lines) embedded in a torus (grey square). The symmetry generators are represented as 
blue circle at the origin (0, 0) for the 7r-rotation generator r, and as a blue dashed arrow - glide generator 
g. Red and green dots represent the identified singular points of the orbifold. The green and red arrows 
indicate the identification of the edges. The centre of the rectangle (black dot) represents the Pkk parity 
centre. On the left panel we show the symmetric case £ = 1 (R4 = -R5), on the right panel the asymmetric 
orbifold with £ = 4/3 (3i? 4 = 4R 5 ). 



2.1 Definition of the "flat" RP 2 orbifold 

The flat real projective plane is described as a quotient space of a plane M 2 modulo a discrete 
symmetry group^j 

T Rp2 = (r,g\r 2 = (g 2 *r) 2 = l) (2.1) 

where r is the rotation of 7r-degrees around the origin of the co-ordinate system, and g is a 
glide-reflection: 

r . f z/4 ~ r(jte) = -va . f y± ~ g(yd = ~Vi + nR A , 2 2 , 

' 1 2/5 ~ r(jte) = -y s ' ' \ 2/5 ~ g(Vb) = Vs + tt^5 

Translations along the two directions are generated by the glide, therefore the space can 
be seen as a subspace of a 2-dimensional torus with radii R4 and R5. In this sense this 
compact space can be thought of as a patch that covers the entire torus once replicated by 
the symmetries in T RP 2. Fermions can therefore be defined here in the same way as they are 
defined on a torus. One can imagine the fundamental domain of the real projective plane, 
obtained by applying the above identifications of points of 1R 2 , as a rectangle whose sides are 
7ri? 4 and ttR 5 long, and where facing sides are identified in opposite directions, see Figure [Tj 
It is easy to see that this orbifold has no fixed points nor fixed lines. While the rotation r 
leaves the corners of the fundamental rectangle fixed, opposite corners are identified by the 
glide: (0,0) ~ (vri? 4 , tcR 5 ) and (0, irRs) ~ (ttR^, 0). The two points, that we label and tt 

2 Note that the structure of the group is entirely defined by the relations between the generators. The 
particular representation of the generators, in terms of isometries acting on the plane, is not necessary but 
helps in visualisation. 
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respectively, are two conical singularities with deficit angle tt, where the metric is still finite. 
In general 4-dimensional Lagrangians can be localised on such singular points, and they are 
required by the divergences of loop corrections: however, the identifications provided by 
the glide will make sure that some residual symmetries of the extended Poincare group will 
survive and play the role of KK parity. 

Note that a second glide can be defined as the combination of glide and rotation: 

9' = <? * r : ( Vi ~ = y ' + *\ (2.3) 

I Vh ~ g{ys) = -2/5 + ttR 5 k j 

We can therefore equivalently define the space in terms of the two glides g and g' . 

The space has two free parameters, i.e. the two radii R4 and R§. In the following we will 
use a convenient parametrisation: calling i? 4 always the largest one of the radii (R4 > R5), 
we define 

m KK = ±- = \ and £= ^ > 1 . (2.4) 

In this parametrisation, ttikk w iH set the mass scale for the KK resonances, while £ will 
parametrise the asymmetry of the rectangle, with £ = 1 describing the symmetric case 
(square). 



2.2 Symmetries of the "flat" RP 2 

Geometrical symmetries of the fundamental space of the orbifold may be interpreted as 
parities or global symmetries acting on the Kaluza-Klein modes of the model. However, 
this is not always the case, and one needs to be careful with the relation between such 
symmetries and the symmetries defining the orbifold. In what follows we will keep denoting 
the generators of the symmetries defining the orbifold by lowercase letters r and g (and </), 
while we will use capital letters for the symmetries of the fundamental domain. 

The fundamental domain of the RP 2 is a rectangle with sides irR^ and irR$. The corners 
of the rectangle are special points corresponding to conical singularities of the space where 
localised counter-terms emerge: in the fundamental space there are only two such points as 
facing corners are identified by the glide and therefore are identical points. Also, opposing 
sides of the rectangle are identified, as shown in Figure [TJ 

The fundamental space is invariant under a rotation by 7r degrees around the centre of 
the rectangle: this symmetry maps the rectangle onto itself, including singular points and 
sides. In terms of the co-ordinates: 

Rf . f J/4 -»■ R'{Va) = -2/4 + 71" #4 / 2 5 x 

' \ y 5 -»■ R'(y 5 ) = -y 5 + 7rR 5 ■ 1 ' ' 

As this symmetry changes sign to both co-ordinates, the components of the 6-dimensional 
fermions with opposite chiralities will necessary pick up opposite parity: therefore, this 
symmetry is not a convenient definition of a parity, because it does not allow do define a 
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single parity to the massive KK fermions. However, one can define an equivalent symmetry 
by combining it with the rotation in the orbifold symmetry group: 



P KK = R'*r:\y^ P p KK ^\ = y ^% . (2.6) 

Now, the symmetry is a simple translation along both directions, and it can be easily seen 
that all KK modes in the tier (k,l) will pick up the same phase (— l) k+l . As such, this is a 
perfect example of KK parity: states with k + 1 odd will not be able to decay into SM states 
and the lightest one, belonging to the (1, 0) tier will be exactly stable^} The symmetry does 
not rely on any assumption on the Lagrangian terms localised on the two singular points, 
therefore it is also respected by the UV completion of the model. 

Imposing conditions on the localised Lagrangians may enhance the symmetries of the 
model. For instance we may assume that the Lagrangians localised on the two singular 
points are the same: this is true, for instance, for the counter-terms generated by loops 
involving bulk couplings. In this case, the space is also invariant under a mirror symmetry 
with respect to the axes y 4 = 

M 4 :l y ^f;[ y4 \ = - y * + 7TR * . (2.7) 
\ y 5 ->■ M 4 (y 5 ) = 2/5 

Once more, one can remove the minus sign in the transformation of y& by combining it with 
the glide g: 

This is a good parity because all states in the tier (k, I) will pick up the same phase (— 1) . 
If this parity were conserved, then transitions like (0, 1) — > (1, 0) or (1, 1) — > (0, 0) would be 
forbidden. This implies that processes that violate P' KK can not be mediated by bulk loops, 
and only by (unequal) localised Lagrangians. 

The case R4 = R5 (£ = 1, degenerate RP 2 ) also seems to have enhanced symmetry: the 
masses of symmetric tiers (k, I) and (I, k) become equal at tree level and the fundamental 
domain becomes a square. The square is invariant under a mirror symmetry with respect to 
the diagonal: 

M d A y ^^\ = y * . (2.9) 

Such symmetry exchanges symmetric states, therefore the spectrum seems to achieve an 
enhanced symmetry. However this is not the case: to see this, we can notice that the 
exchange of the two directions will also exchange the two glides, g -H- g' . Fields which 
have different parity under the two glides {p g i = p g p r ), therefore, will not be mapped into 
themselves by this symmetry: this always happens for fields odd under the rotation. To 
see this explicitly, let's consider the KK expansion of a scalar field with parities p r = — 1 



3 In the case R4 — R 5 , there may be two degenerate DM candidates, as the (0, 1) is degenerate with (1, 0) 
at three level. More details will be discussed in the following sections. 
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and p g = 1 (therefore, p g > = —p g ). the spectrum contains tiers (2k, 0) but not (0,2k) [2], 
therefore the map cannot be closed. Similarly, tiers (0, 2k — 1) do not have a symmetric 
state. Even if one choses to include only bosons with parity p r = +1 (for which p g > = p g ), 
problems arise when considering fermions. In fact, the left- and right-handed components of 
the 6D fermions do have opposite parities under the rotation, therefore it is unavoidable to 
have fields with p r = — 1. Nevertheless, the spectrum is symmetric and gauge interactions 
do not depend on the p g of the field. On the other hand, Yukawa couplings are proportional 
to p g , therefore the symmetry will flip the sign of the Yukawa couplings. There are other 
localised interactions that may distinguish between the two 6D chiralities and depend on p g . 
For such reasons, cannot be considered as an exact symmetry of the orbifold. 

To summarise our discussion, we showed that, independently on the value of the radii, 
the RP 2 possesses a unique exact KK parity Pkk = (—l) k+l - 

In the following, we will consider the space-time Ai 4 x RP 2 with a 6-dimensional flat 
metric g = diag(l, —1,-1,-1,— 1,-1). We denote by x M the co-ordinates of .M 4 with 
/i G {0, 1, 2, 3} and by y a the co-ordinates on the orbifold with a G {4, 5}. When needed we 
note a vector y = (1/4,1/5). The capital letters of the alphabet run over all the co-ordinates 
of the six-dimensional space, M, N G {0, 1, 2, 3, 4, 5}. 



2.3 Parities of the fields and KK decomposition 

In this section we summarise the Kaluza-Klein decomposition of quantum fields defined on 
the six- dimensional background A4 4 x RP 2 described in the previous section. We introduce 
the fields defined in six dimensions: a scalar field <&(x^,y a ), a gauge field A M (x fJ, ,y a ) and a 
spinor ^>(x^,y a ). For every orbifold symmetry, namely the rotation r and the glide g in our 
case, it corresponds a transformation of the fields: 

X(r(y a )) = Pr nX(y a ) and X(g(y a )) p g QX(y a ) , (2.10) 

where 1Z and Q are the transformation operator in the adequate spin representation and 
X stands for any quantum field. The invariance under the orbifold projection requires to 
impose such properties on each field in the theory, where the two parities (p r ,p g ) determine 
the spectrum and property of the field. From the properties of the symmetries r and g in 
Eq. |2.1 , it can easily be shown that the possible parities of the fields under the rotation and 



glide-reflection are (p r ,p g ) = (±, ±) respectively^} Thus every field will be characterised by 
a pair (p r ,Pg)- As the RP 2 has no boundaries, we can decompose each field in KK modes as 
on the torus, where the wave functions are combinations of sine and cosine functions of the 
extra co-ordinates, and then impose the orbifold projection on each field. 



4 The glide can also pick up an imaginary phase, thus (p r ,Pg) = (ij i*) are also allowed choices. We will 
not consider them here for simplicity. 
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Scalar field 



As a simple example, we can detail the case of a scalar field, for which 1Z and Q are equal 
to the identity operator. The 6D action of a complex scalar field &(x fl ,y a ) of mass M<j, isQ 

dV / dy A / dy 5 [D M &D M $ - M%&$] (2.11) 
Jo Jo 

with Dm = 9m — W^liK the covariant derivative. The field is decomposed into KK states 
labelled by a pair of positive numbers (k, I), with k, I > 

*(*",!/") = X> (W 0O/w(jO, (2-12) 

fc,Z>0 

where the wave function is generically a combination of sin and cos: 

fk,i(v a ) K [sinfe/^) ,cos(ky4/R4)) x [sin^/i^) , cos(Z2/ 4 /i? 4 )] . (2.13) 

Solving the equations of motion, one finds the mass spectrum of the Kaluza-Klein resonances 
at tree level : 

{-d\ - dl + Ml) f k>l (y a ) = My kjl {y a ) = (m^ + M|) f kil (y a ) (2.14) 

with 

m h = %> + 4 = V^ 2 + e/ 2 . (2.15) 

The functions f k ,i are found by imposing the parities (p r ,p g ) on the wave functions, so 



5 In this action we integrate over the fundamental domain of the torus. Integrating over the domain of 
the RP 2 would only change the normalisation of the fields, however the physical properties, and the KK 
expanded Lagrangians, are exactly the same. 
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there are four possible KK towers with different parities: 

oo 

= + E (cos(2% 4 / J R 4 ) 0( 2fc <°) + cos(2% 5 / J R 5 ) 0(°' 2fc )) + (2.16) 

fc=i 

+ ^ ( COS (%4/^4) COs(ly 5 /R 5 ) 4+?=2m + 8i n (W#4) SUl^/i^) 0i+?=2m+l) > 

fc,«>0 

$ (+ - ) = ^ ^ (cos((2fc - l)y 4 /i? 4 ) 0^^^ (2.17) 
fc=i 

+ ^ ^sin(% 4 / J R 4 ) sm(ly 5 /R 5 ) 4>t+i=2 m + cos(% 4 /i? 4 ) cos(ly 5 /R 5 ) ^+f =2m+1 ) , 

fc,Z>0 

oo 

$ ( - +) = V (sin(2% 4 /i? 4 ) (2fe ' o) + sin((2A; - l)y 5 /i2 5 ) (O ' 2fe_1) ) + (2.18) 

v 2 ^ fc= i 

+ ]V S ( Sin (%4/i? 4 ) COs(ly 5 /R 5 ) 4> ( k+l=2m + COs(% 4 /i? 4 ) Sm(ly 5 /R 5 ) 0fe+?=2m+l) > 
fc,«>0 
1 oo 

$ ( " } = ^ V (sin((2A; - l)y 4 /i? 4 ) (2.19) 
^ 2iV fc =i 

+ Jf Yl ( cos (%4/^4) sin(/y5/ J R 5 ) ^ =2m + sm(ky 4 /R 4: ) cos(ly 5 /R 5 ) 0i+?=2m+i) 5 

fc,«>0 

where TV = n^/R^R^. The only case with a zero mode is (+, +), therefore a model with bulk 
Higgs must have such a parity in order to include a Higgs boson in the zero mode spectrum. 

Gauge fields 

Gauge fields in 6 dimensions are given by a 6-component vector: 

A M (x^y a ) = (A^,y a ),A^,y a )). 

Being a vector, the parities of the various component of Aj^ix^, y a ) are sensitive to the 
action of the symmetry on the co-ordinates. In fact, under rotation and glide, a 6- vector 
transforms as 

A M (r(y a )) = K(A^y a ), -A 4 (y a ), -A 5 (y a )) , A M (g(y a )) = G(A^y a ), -A 4 (y a ), A 5 (y a )) , 

where TZ and Q are suitable gauge transformations. In the following we will focus on a 
model where the gauge groups are the same as in the SM, therefore we will require that all 
the gauge fields have a vector zero mode, associated with an unbroken 4D gauge symmetry. 
Therefore, for all the gauge groups, it is enough to chose 1Z and Q equal to the identity. More 
in general, a non trivial gauge transformation will allow to break the 6D gauge symmetry to 
a subgroup in 4D (namely, not all the gauge bosons will have a massless vector mode, but 
only the gauge fields associated with a subgroup). 
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For a (+, +) vector, which will contain a zero mode vector boson in the KK theory, the 
parity of the "scalar" components will be (— ,— ) for At and (— ,+) for A 5 . The spectrum 
of gauge fields is complicated by the mixing between the vector components and the scalar 
ones, and an arbitrary gauge choice. A complete and detailed analysis of the spectrum can 
be found in [2]. Here we will limit ourselves to a brief description of the spectrum in the 
Feynman-'t Hooft gauge, where the goldstone bosons appear explicitly in the spectrum with 
the same mass as the vectors they belong to. The spectrum of the A^, v4 4 and A 5 components 
is the same as for a scalar field listed before. The massive vectors in the spectrum pick up 
their longitudinal degree of freedom by eating up a linear combination of the two scalar 
components. While the zero mode is left massless, the (2k, 0) vector modes eat up the 
(2k, 0) modes in A^ while the (0, 2k) ones eat up the (0,2/c) modes in A4. On the other 
hand, the (2k — 1, 0) modes in A4 and the (0, 2k — 1) modes in A§ are physical scalars 
because there is no corresponding vector boson with the same mass. For the (k, I) modes, a 
linear combination of A 4 and A 5 is eaten, while another one is left in the spectrum. 

Spinor field 

The spectrum for fermions is also complicated by the action of the orbifold symmetries 
on the components of the fermion field, and we invite the reader to consult Ref. [2j for a 
detailed description. Here we mention only some points that are important for our purposes. 
In 6 dimensions it is possible to define two independent 6D chirality operators (which are 
different from the 4-dimensional chirality!). A generic 6D fermion has four components that 
are Weyl spinors: \1/ 6D = (x+, fj-, X-i V+) T an d in this notation ± subscripts correspond to 
the 6D chiralites while Xi V are the 4D chirality eigenstates. The glide symmetry flips the 
6D chirality of the fermion it is acting on as 

p g G^(y a ) = P g G(x+, x-, v+) T = Pg(X-, V+, X+, V-) T , ( 2 - 2 0) 

therefore a consistent model with fermions must contain 6D vector-like fermions. The parity 
under the glide p g does not play any significant role and it's only constrained by the presence 
of Yukawa couplings because the p g of the Higgs field is fixed to be equal to +1. On the other 
hand, a bulk mass, that would connect the two 6D chiralities, is forbidden by the rotation 
symmetry. For a vector-like fermion in 6D, the parity under rotation selects the chirality of 
the zero mode|^j p r = +1 for a left-handed zero mode and p r = — 1 for a right-handed one. 
For each choice of the parities, the KK tower of a fermion will contain a chiral zero mode, 
one vector-like fermion in each tier (k, 0) and (0, k) and two vector-like fermions in each 
(k,l) tier. Note also that for each SM fermion, the 6D theory will contain an independent 
vector-like 6D fermion for the two chiralities: for instance, the electron will be associates 
with a 6-dimensional spinor e| D with p r = +1 corresponding to the left handed component 
of the electron, and a e^P with p r = — 1 corresponding to the right handed component. 

6 Note that the phase under rotation of the left-handed and right-handed components of the 6D fields are 
opposite: p r TZ^(y a ) = p r lZ(x+, X-> V+) T — Pr(x+> ~V-> X-> ~V+) T ■ What we define as the parity of 
the 6D field is the parity of its left-handed components. 
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fl, 1,0) 


( + ,+) 


V 




V 


V 
V 


L 


(1,2,-1/2) 


/ J \ 

(+,P l q ) 


y/l.h. 
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VV 


E 


(1,1,-1) 


(-P l q ) 


V r -h- 
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V 


vv 


Q 


(3,2,1/6) 


(+,Pl) 
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V 


VV 


u 


(3,1,2/3) 


(-P q q ) 


V r.h. 


V 


V 


VV 


D 


(3,1,-1/3) 


(~,Pl) 


V r.h. 


V 


V 


VV 


H 


(1,2,1/2) 


(+,+) 






V 


V 



Table 1: UED spectrum on the RP 2 . The subscript (j> indicate the physical gauge-scalar. The parity under 
glide of leptons p l g and quarks jfl g are arbitrary, and due to the fact that the fermions are 6D vector-like, 
either choice will lead to the same physics. 



UED model 

In the following, we will focus on an incarnation of Universal Extra Dimension (UED) models 
on the RP 2 : in this model, each SM field is the zero mode of a 6-dimensional field with 
suitable parities. The tree-level spectrum of the model is summarised in Table [TJ In the 
degenerate case £ = 1, when i? 4 = R 5 = R, the symmetric levels (k,l) and (I, k) are 
exactly degenerate at tree level. When £ > 1, the degeneracy is removed; moreover, as £ 
increases, there are crossovers between different levels. This fact can be very important for 
the phenomenology of the model both at the LHC and regarding the Dark Matter abundance 
calculations. As shown in Figure [2j the crossover can take place for small values of £ among 
more massive levels. On the other hand, the crossover for light modes requires £ — 1 to be 
order 1: for instance, the level (0, 1) becomes degenerate with (2, 0) for £ = 2 and level (1, 1) 
with (2,0) for £ = ~ 1.7. 

3 Loop corrections to the level (2,0) and (0,2) masses 

Loop corrections play a crucial role for the phenomenology of UED models, because the 
states in each tier are exactly degenerate at tree level (ignoring effects from the electroweak 
symmetry breaking). Furthermore, for even tiers, loops can also induce vertices that mediate 
decays directly into a pair of SM particles, thus violating the KK momentum conservation 
rule. For the lightest odd tier, the mass splitting will determine the typical visible energy 
of the decay chain, which will inevitably end into the Dark Matter candidate; for even tiers, 
the loop corrections are crucial to determine the branching ratios into resonant SM pairs, 
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Figure 2: Crossover between tiers as a function of £. The colour-code indicates the KK parity of the tiers: 
red for odd tiers, blue and green for even ones where the green ones can only decay into zero modes via the 
localised interactions. 



the relevance of Drell-Yan resonant production via a heavy state, which can impact both 
the LHC search strategy and add resonant contributions to the Dark Matter annihilation 
cross sections. The latter determines the relic abundance of Dark Matter and therefore the 



preferred mass range to explain current cosmological observations 10 . 

In the rest of this section, we will introduce the loop corrections to the masses of the 
even states (2,0) and (0,2), focusing on the dependence on the two radii (R4 and R$, or 
equivalently ttikk and £) and the relevance of the mixing between the two tiers. We will also 
recall the analogous formulas for the odd states (1,0) and (0,1) [2], including the dependence 
on the two radii. Loop induced KK-number violating vertices and mixing terms between 
the two even tiers, which also violate KK-number conservation, only depend on the localised 
divergences, thus they can be extracted by using the localised counter-terms as in [5]. In the 
following, we will only calculate the mass mixing in detail. 

Mixed propagator formalism 

The loop corrections to modes with a zero KK number, like (n, 0) and (0, n) can be conve- 
niently computed by use of mixed propagators in 6 dimensions: we will Fourier-transform the 
non-compact co-ordinates into momentum space p M , which contains the physical momen- 
tum and energy measured by 4-dimensional observers, and leave the other two co-ordinates 
in position space. The 6D propagator can in general be expressed in terms of Bessel func- 
tions 11 , which are not very handy in loop calculations. Using the presence of a zero KK 
number for the external lines of the loops, we can expand in KK modes along the direction 
corresponding to the zero external KK momentum. For instance, for modes (n,0), we can 
expand in KK modes along the second co-ordinate y 5 . The 6D propagators can therefore 
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be written in terms of 5D propagators: for instance, for a scalar propagating on a torus, we 
have 

oo 

GT S (>, y- y>) = J2 G r lc (», |V4 - v'A) fito)M) • (3.21) 

l=— oo 

In this expression, /;(y) are the normalised wave-functions on the circle and / labels the 
masses of KK states expanded along y 5 , 

M =vm- (3 - 22) 

and G$ rcle are the 5D propagators jl2| of a scalar field of mass 

^drcic/ i m • cos xii^Ri - |s/4 - y'A) , , ■ ,,, 

(»,ltt-ttl)=» 2x;Sin(x ^ 4) » Wlth K = \h 2 -^- ( 3 - 23 ) 

The advantage of this expressions is that the integrals on the co-ordinate y 5 are very simple 
as the wave functions of the external states do not depend on it, thus we are left with 
integrals of exponential functions. The propagators of vectors and fermions can be written 
in terms of the scalar propagator: in Feynman-'t Hooft gauge 

GaT = -J M C , and G^ us = (pT^ - iT 4 d y4 - iT 5 d ys ) G$ orus , (3.24) 

where are the Dirac Gamma matrices in 6D. 

Finally, the propagator on the orbifold can be written in terms of the torus propagator 
observing that the propagators from y to the images of y' via the symmetries of the orbifold 




are physically identical. Therefore, one obtains 13 



G% p2 (f, y, y') = \ [GT S (>, y - y) + v 9 GT S (p",y- g{rh) 

+ P 9 Pr GT S {f,V-\9* r] (y')) + Pr GT ls (f, y - r (y')) ] ,(3.25) 

where {p g ,p r ) are the parities of the propagating 6D scalar field. The factor of 1/4 takes 
into account the fact that we are integrating over the fundamental space of the torus, which 
is 4 times bigger than that of the RP 2 . Note that all the wave functions are also normalised 
over the torus, without affecting the physical results. For vectors and fermions, we need to 
additionally multiply by the transformation operators in the appropriate spin representation, 
thus, for instance 

/ . „MN r / \ / \ 

gZ 2 = -— [GT S (f,y-v') + (-i)^P 9 G<r us \fj-g{$)) 

+ (-1)^)^^^^-^,.^) +(-l)^)+V.5 )pr G' t - us ^,y-r(|?)) • 

(3.26) 
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General loop structure 



Using this form of propagator as in Eq.(3.25), a generic loop correction to any field can be 
decomposed as 

iU(m KK , = 1 (~n.T(m KK , f ) + p g U G (m KK , f ) + p g p r Yi GR {m KK , f ) + p r U R (m KK ,^)) , 

(3.27) 

where the parities refer to any one field propagating in the loop. 



The first term in Eq.(3.27), Ut, corresponds to the loop corrections of a theory defined 
on a torus: this contribution is generically divergent, however it gives a finite contribution 
after renormalisation of the bulk kinetic terms. Here we will follow the prescription of 



Ref.s 14, 15 , where we remove the contribution of the zero winding modes, i.e. the modes 
that do not wrap around the torus. The terms Uq and Hqr correspond to glide symmetries 
and are finite. The reason behind the finiteness is that the points y and its image g(y) (or 
[g * r](y)) never coincide, i.e. the glide(s) do not admit any fixed points. The last term, ITr, 
corresponds to the rotation transformation, and it is divergent because the rotation admits 
fixed points. It was shown explicitly in [2] that divergences only arise in the points where 
y* = r {yl), i-e. on the corners of the rectangle. Thus, such divergences can be renormalised 
by adding counter-terms on the two singular points, whose structure has been studied in (5] . 
The bulk loops simply require equal terms on the singular points, as loop interactions cannot 
distinguish between the two. 



The loop contributions generally depends on both radii, i.e. ttikk and £ as in the Eq ]3.27 
One can easily relate the loop contributions to the modes (n, 0) and (0, n) by simply exchang- 
ing the two radii, and the two glides: 

p g p r U^ \Cm K KA/0+Pr^ 0) ^m KK ,l/0) . (3.28) 



From Eq. 3.28 it is clear that, in the degenerate case £ = 1, we have n>>°) = IT '™) only if 
Pa = Pg' = VgVr (i-e. p r = 1) for all the fields running in the loop because, as we will see, 

iWn GiJ . 

The loops will also generate mass mixing between states with different KK numbers, as 
long as they are allowed by the symmetries of the space. For instance, states (rii,0) and 
(n 2 ,0) can mix via a one loop diagram of modes ( ni + n2 , m) and ^ l Wl ~" 2 l ; only if both 

ri\ and n 2 are even or odd. Analogously, modes (rii,0) and (0, n 2 ) can mix only for even n\ 
and 7i2 (which are allowed to mix with n = 0). The mixing effect is negligibly small, unless 
the modes that mix are degenerate at tree level: in the following we will be particularly 
interested in the mixing between (2, 0) and (0, 2) which are degenerate in the equal radii 
£ = 1 case. The odd modes, like (1, 0) and (0, 1) cannot mix via loops, however they can mix 



via asymmetric localised terms. From Eq. (3.27), we can infer that the off-diagonal mixing 



terms can only be generated by the rotation-term H R , because both the torus and the glides 
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preserve KK number [J] 



3.1 Gauge bosons 

In this section we will present the results for the mass corrections to the tiers (n, 0), including 
novel results for even n (vectors) and summarising the results for n odd (scalar) from (2j. 



The results for (0, n) modes can be extracted by use of Eq. 3.28 The loop corrections to the 
masses can be written in general as: 



a™ 2 - g2C ( r ^ i A 

dm (n,0) - 64?r 4 m KK A (n,0) , 



(3.29) 



where g is the gauge coupling, C(r) the Casimir of the representation r of the field running 
in the loop (to be substituted by the charge squared for a U(l) gauge boson), and the loop 
coefficients Ar n ^ are given in the following table: 



An,0) 


torus 


glide (p g x) 


glide' (pgPrX) 


rotation (p r x) 


, n-even 
gauge loops , , 
n-odd 


4T 6 (0 


14C(3) 


e(U((3) + V g (n/0) 
e (14C(3) + S g (n/0) 


8n 2 n 2 L 


r . , n-even 
fermion loops , , 
n-odd 


-8T 6 (0 











, , n-even 
scalar loops , , 
n-odd 


mo 


7C(3) 


e 2 (7C(3)+y s (n/e)) 

e (7C(3) + S s (n/0) 


-|n 2 7r 2 L 
n 2 n 2 L 



where C(3) is the Riemann Zeta function, 



1 



E 



1 



7r^ ^— ' (n 2 + m 2 /^ 2 ) 2 

{n, m} # {0, 0} 



(3.30) 



encodes the contribution of the torus loops (numerically, T 6 (l) = 1.92. ..), L — log(A 2 /m 2 KK ) 
encodes the logarithmic dependence on the cut-off A, and the functions S g>s and V 9jS contain 
small and finite corrections due to the glides. The latter functions and T 6 (£) can be found 
in the Appendix |XJ From the table it can be seen that the difference between the scalar 
(n-odd) and vector (n-even) components of the gauge fields comes from the contribution of 
small finite glide terms, and of the logarithmic divergent contribution of scalar loops. 

In the UED model we are interested in here, the only bulk fields are the SM ones, 
with parities (+, +) (except for fermions, whose parities do not enter the mass corrections), 
therefore one can calculate the mass corrections for the U(l), SU(2) and SU(3) gauge bosons 

7 In fact, the orbifold defined by the glide alone is the Klein bottle, whose loop contribution is given by 
IlKicin — 5 (IIt +p 9 IIg). The Klein bottle is invariant under translations along the two directions and has 
no fixed or singular points, therefore no KK-number violating mixing can be generated by either IIt or IIg 
(and U GR ). 
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as follows: 



^rriKK [-36T 6 (0 + (1 + £ 2 )42C(3) + 3£ 2 S 9 (n/0 + 24n 2 vr 2 L] 
^i m KK [-79T 6 (0 + (1 + e 2 )7C(3) + e 2 5 s (n/0 + nVL] , (3.31) 
gfljmmr [-39T 6 (0 + (1 + £ 2 )35C(3) + £ 2 (2S 9 (n/£) + S s (n/0) + 17n 2 vr 2 L] , 
for the scalars ra-odd; and 



S< 4 = # " a 
5m k = 



5mL = 



= 6^4^[-36T 6 (0 + (l + e 2 )42C(3) + 3^V s 

_2 



r 2 Pi 2 

64vT 4 

r 2 ^?2 2 

5m w = -zr—^m, 



n/i) + 24nVL] 
79T 6 (0 + (1 + ^ 2 )7C(3) + CV a (n/0 - ^nVL , 



64vr 4 m ^ 



(3.32 

O J 

47 

-39T 6 (£) + (1 + e 2 )35C(3) + e(2V g (n/£) + F s (n/£)) + y™ 2 ^ 



two 



for the vectors n-even. In the n-even case, loops will also induce a mixing between the 
tiers (n, 0) and (0, n) proportional only to the divergent rotation contribution. Such terms 
can be computed using the localised counter-terms, which are kinetic terms for the vector 
bosons [5]. The off-diagonal contribution to the mass can be written as: 



r 2 9 2 C(r) 2 1 + £ L 

5m (n,0)-(0,n) = 64?r 4 m ^ A M) ' ( 3 -33) 

where the A< n ,a) coefficients are the same as in the table (for the rotation only). For the 
UED model, the off-diagonal terms are 



where the At 



5m 



5m 2 , 



= A_m> i±i! r 2 4n 2 7r 2 Ll , 
(n,0)-(0,n) 64vr 4 KK 2 1 J ' 

^2 1 , t 2 r 1 i 



5m 



U' u 



<?? 2 1 + e 2 

— ± 777, 

(n,0)-(0,n) 64tt 4 ^ x 2 

9l 2 1+e 2 

— _ 771 

(n,0)-(0,n) 647T 4 KK 2 



1 2 2 r 

— n 7r L 
3 

'47 2 2r ' 
— n 7r L 
3 



(3.34) 



£ dependence 

In Figure [3j we plot the relative value of the loop corrections to the 3 gauge bosons of U(l) 
(B), SU(2) (W) and SU(3) (G) for the tiers (1, 0), (0, 1), (2, 0) and (0, 2): the plots only show 
the diagonal loop corrections, so the effect of cross-level mixing and of the Higgs VEV are 
not included. The three colours correspond to different values of the cut-off A, in particular 
in green we show the correction for AR = 1, i.e. for vanishing divergent log (L = 0). For 
the odd levels (on the left), the plots show a strong dependence on £ for the lightest mode 
(1,0), while the dependence is much milder for the (0,1) level. The reason behind this 
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Figure 3: 8m 2 g /m 2 g as a function of £: the three colours correspond to Ai? = 10 (L = 4.6), Ai? = 5 (L = 3.2) 
and Ai? = 1 (L = 0); the solid line to tier (n, 0) while the dotted one to tier (0, n); in the first column n = 1, 
in the second n = 2. 

behaviour is to be traced back to the strong £ dependence of T 6 (see the Appendix [A] for 
more details) that contributes to the finite part of the correction, while the divergent part 
is independent on £. In fact, for £ > 1 the correction tends to more negative, due to the 
negative contribution of the T$ term. In the (0, 1) case, dividing by the KK mass ^ 2 m 2 KK 
cancels out the growth with £ of the finite part, thus the result is almost ^-independent. The 
finite part always plays a crucial role for the hypercharge, because of the smallness of the 
divergent contribution which is coming only from the Higgs loop. For the even levels, this 
feature is even more obvious: the loop corrections are dominated by the divergent parts, 
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thus the ^-dependence is suppressed, except for the U(l) gauge vector B. Nevertheless, the 
£ dependence is always milder than for the odd state. This feature repeats itself for modes 
with larger n, both even and odd: the point is that the divergent term is proportional to 
m 2 x , while the finite contribution are only proportional to 1/R 2 and therefore are suppressed 
by a factor 1/n 2 with respect to the divergent ones. 



Tier mixing 

For the odd modes, no mixing takes place, unless non negligible and asymmetric localised 
terms are present, therefore masses are well defined by equations |3.32 The localised terms 



always correspond to higher order operators (dimension 8 compared to the standard dimen- 
sion 6 bulk terms) (5), therefore we may expect them to be small compared to loop induced 
corrections. 

In the (2,0)- (0,2) system, the mixings between the levels cannot be neglected because 
the off-diagonal terms are generated at loop level. For simplicity, we will define the loop 
corrections to a generic A( n ^ and A(o, n ) gauge vector (with n even) as 

1 + £2 

6m A\(n,0) = n2m KK^(n,0) , 5m A\ { ^ n) = n2 ^ m KK^(0,n) , 5m 2 A \ (n)0) _ (M = ^ '^—m 2 KK 6' n . 

The two new mass eigenstates are defined as 

A;- = cos <^4 n ' 0) -sin y n Af n \ (3.35) 
A; + = sincpnA^ + coscpnAf^; (3.36) 



with masses 



where 



m\ ± = n 2 m 2 KK 



(l + e 2 + 5 (n ,o)+^(0,n) ±A n) ( 3 - 37 ) 



a„ = ^(i + e)w + (e - 1 + - 5 (nfi) ) 2 (3.38) 

and the mixing angles (p n are a non-trivial functions of £: 

tan Vn = ^-(f-l+^W-W . p.39) 
Note that for £ 2 — 1 3> 5's we have 



m 2 n± = n 2 m 2 KK x I / (3-40) 



with eigenstates A n ~ ~ A^' ^ and A™ + ~ A^' n \ therefore the mixing will be negligible and 
we can and the two mass eigenvalues will be proportional to itlkk and ^tukk- 



18 



In the degenerate case £ = 1, the loop corrections to the two modes are the same, 
<5(n,o) = <5(o,n) — &n'- this is true because the UED model does not contain bosons with 
p r = — 1 (let us stress again here that this is not a generic feature of the RP 2 ). In this case, 
A n = 25' n , and ip n = it /2: the mixing is maximal and the mass eigenstates are given by 

< = ^(4"* 0)± 4°' n) ) ( 3 - 41 ) 

with masses 

m 2 n± = n 2 m 2 KK (l + 5 n ±5' n ). (3.42) 

As the divergent part in 5 n and 5' n have the same coefficient, it turns out that the fields A™~ 
are insensitive to divergences: this can also be explained by the fact that their wave function 
vanishes on the singular points. 



3.2 Fermions 



In the fermion case, we can express the one-loop corrections from gauge loops to the masses 
as: 



64vr 4 



-m~ KK F(„ j0 ) , 



(3.43) 



where the fermion is in the representation r/ and the group theory factor 62(77) is replaced 
by the charge squared for U(l) loops. For scalar loops via Yukawa interactions, it is enough 
to replace fi^C^r^) — > y 2 , where yj is the effective 4D Yukawa coupling. The coefficients 
-F( n ) are summarised in the following table: 



F(n,0) 


torus 


glide (p g x) 


glide' (p g p r x) 


rotation (p r x) 


, n-even 
gauge loops , , 
ra-odd 







14C(3) 


e(7((3) + F a (n/0) 


An 2 n 2 L 


, . n-even 
scalar loops , , 
ra-odd 







7C(3) 


em3)+F s (n/0) 


\n 2 -K 2 L 



where the small and finite functions F 9tS are shown in Appendix [A], and the parities refer to 
the parity of the bosons in the loop (gauge or Higgs). The only difference between even and 
odd n is in the contribution of the glide, which is absent for even tiers. 

In the UED model all parities are (+, +), thus the mass corrections for a generic fermion 

is 



5m 



M)ln- 



mKK 



E ^(r/) j (e 2 7C(3) +eF g {n/i) +An\ 2 L) 
1 



\gauge 



+y 2 V?C(3) + e F s {n/0 + -nVl 



(3.44) 



5m 



(n,0) 



n—odd 



m KK 
647r 4 n 



E 9iC 2 (r f ) ((2 + £ 2 )7C(3) + eF 9 (n/0 + An 2 rr 2 L) 



\gauge 
2 



2/2 ^7C(3)+£ 2 F>/0 + ^VL 



(3.45) 
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Figure 4: Smf/irif as a function of £ for singlet leptons: the three colours correspond to AR = 10 (L = 4.6), 
Ai? = 5 (L = 3.2) and AR =1 (L = 0); the solid line to tier (n, 0) while the dotted one to tier (0, n); in the 
first column n = 1, in the second n = 2. The plots show a mild dependence on £ in both cases. 

The n-even fermions mix via localised counter-terms: in the fermionic case, there are both 
kinetic terms and bulk mass terms (i.e. terms with a derivative along the extra co-ordinates), 
however the off-diagonal mass corrections can be cast in a simple way 

3 m (n,0)-{0,n) = ^ Q^i^ m KK F{nfl) , (3.46) 
5m(p,n)-(n,0) = £ Sm( nfi )-(0,n) , (3.47) 

where only the divergent terms in -F( n ,o) contribute. Note that the term proportional to £ 
may be <5^(n,o)-(o» ; depending on the parity of the field: in fact, the off diagonal mass terms 
are proportional to the mass of the fermion component which is odd under the rotation r, 
the reason being that the non vanishing localised term must be proportional to the derivative 
of its wave function. Therefore, with the standard mass definition in 4 D mipip, where ip is 
a generic Dirac spinor describing the fermion KK modes, then 5m.(o in )-(n,o) is proportional 
to £ for a fermion with left-handed zero mode, while the other situation occurs for fermions 
with right-handed zero modes. 
For a UED fermion: 

rriKK 

5m ( n,0)-(0,n) = 

£ dependence 



4 E 9lC 2 {r f ) + \y) 

gauge 



n 2 7r 2 L. (3.48) 



In Figure |4J we show the relative mass correction 5m /m for the fermions. In the plot we 
show the results for the singlet leptons E, as the behaviour is very similar for other fermions. 
The plot shows a very mild dependence on £ for both even and odd tiers: this behaviour is 
to be understood in relation to the absence of bulk contributions to the mass proportional to 
T 6 (£). The masses, therefore, are always dominated by the divergent terms, while the finite 
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contributions are always very small (green lines in the plot). The plot also shows that the 
relative mass correction is very similar for the two modes (n, 0) and (0,n). 



Tier mixing 

Similarly to the bosonic case, the only mixing generated at loop level occurs for the even 
modes. For a generic value of £, one can diagonalise the mass matrix by use of two rotation 
matrices, acting on the two chiralities of the fermion (labelled by L and R for left- and 
right-handed chirality respectively): 

1pn-,L/R = COS ip n ,L/R 4>(n,0),L/R ~ SU1 <p n , L / R ljj( ,n),L/R , (3.49) 
^n+,L/R = Sin <£ n)L / R 4>(n,0),L/ R + COS <p n , L /R 1p(0,n),L/R ■ (3.50) 

For a KK tower with a left-handed zero mode, the mass corrections can be written as 

^ m V>l(n,0) = nm KK(l + S n ) , <fotyl(0,n) = n i m KK^\ + S n ) , 

<WI(n,o)-(o,n) = nm KK 8' n , <Jm^| (n>0) _ (Oin) = nim KK 5' n ; 

where we have assumed that the loop corrections for the two tiers are the same, $< n ,o) ~ 
<^(o,n) — &n- I n this approximation, the mixing angles tp n ,L/R and the mass eigenstates m| ± 
and are given by 

ta "^ = 2(1 + 4.) • (3 ^ 51) 

tan ^« = 4^(1 + *,) ; (3 ^ 52) 



and 



where 



2 n 2 l 



t = ^2 ( (1 + ^ ( (1 + 6n)2 + ^) ± A «) ' (3 - 53) 



= V(! + £ 2 ) 2 ((1 + V + « - 4^ ((1 + 5 n y - <^ . (3.54) 

As before, for £ — 1 S' n , the two tiers decouple. For a KK tower with right-handed zero 
mode, the only difference is that 

6m ^\(n,0)-(0,n) = < m KK5' n , <frjty | (n>0) _ (0>n ) = 5 

and the two mixing angles are exchanged, while the mass eigenvalues are unaffected. 

In the degenerate radii case, £ = 1, the situation is similar to the bosonic case: the mass 
corrections for the two degenerate tiers are exactly the same, and the mixing between mass 
eigenstates becomes maximal. In fact, for £ = 1, A^ = 4(1 + 5 n )5' n , and = (f n ,R = 7r/2. 
The mass eigenstates are the sum and difference of the two modes, with masses 

m n ± = nm KK (l + 5 n ± S' n ) . (3.55) 

As in the bosonic case, the mode ip n - only received finite corrections due to the fact that 
the divergent part in 5 n and 5' n are the same. 
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3.3 Scalars (Higgs field) 

In UED models the Brout-Englert-Higgs mechanism is implemented in the same way as in 
the Standard Model, i.e. by introducing a bulk scalar doublet which, due to a properly 
engineered potential, acquires a vacuum expectation value (VEV) and gives mass to W, Z 
and fermions. Here we will focus on the simplest possibility that the Higgs has a zero mode, 
thus parities (+, +), and that the VEV is generated by a negative bulk mass. Denoting by 
&h the bulk scalar, the tree level potential is 

\ 2 

£buik = {D^ H ) ] D^ h - m 2 B & H $ H + -f ($ f H $ H ) , (3.56) 

where me is the bulk mass, As is the quartic coupling; note that the quartic operator has 
dimension in mass 8 (as the scalar field has dimension 2), therefore the coupling A# has 
dimension —2. This fact may rise the problem that, in order to have a heavy enough Higgs 
boson, we need a higher order operator with unnaturally large coupling: here we will brush 
this issue aside, as the Higgs sector of this model should just be taken as a toy model because 
it has more serious issues than the SM Higgs itself. 

Loop corrections will also affect such potential and their structure is similar to the gauge 
bosons: there are finite contributions from the torus and glides which are almost mode- 
independent, the only difference coming from the small n-dependent contributions generated 
by the glide', while the log divergent terms are proportional to the mass of the mode. The 
loops also contain quadratically divergent contributions which correspond to a localised 
mass. The mode-independent corrections can be reabsorbed in the bulk mass, which is a 
free parameter, therefore we will neglect them in the following: the only contributions that 
affect the spectrum are thus the localised contributions that can be parametrised in terms 
of the counter-term Lagrangian [5] 

+ c' 5 (dl& H $> H + & H 9 5 2 $ H ) - 5m 2 & H <S> H , (3.57) 

with coefficients 

3 $ + ^ oS ' (3 ' 58) 

Note that 5m 2 , which is a brane localised mass, is a free parameter, and that the radiative 
corrections do not include loops proportional to the quartic coupling. Furthermore, the 
localised mass term is a dimension 6 operator (4 from the two scalar fields and 2 from the 
localising delta functions), therefore it must be considered as a new tree level parameter: 
not that the coefficient 5m 2 /A 2 is dimension-less. 
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In order to understand the spectrum, we need first to study the effective potential for 
the zero mode, which corresponds to the SM Higgs boson. In terms of the bulk and localised 
Lagrangians, the potential for the zero mode is 

V((f>(o,o)) = (m| + m 2 b )4>\ m <p {m + ^A(0{ o o) 0(o,o)) 2 , (3.60) 

where A = 47r2 ^f rfi6 is the effective 4D quartic coupling and ml = 2 -k^r.r 6 ls contribution 
from the localised mass. After the usual expansion 



(3.61; 



we obtain that the SM VEV v and the Higgs mass are given by 

m 2 H = -2(m| + m 2 b ) = Xv 2 . (3.62) 
In the following, we will consider m^ a free parameter and replace mg with m# (a measured 



quantity) in all expressions by use of the relation in Eq. 3.62 Note also that we can expect 
m b to be below the TeV scale in order not to introduce a severe fine tuning in the Higgs mass: 
this is the only concession to naturalness we will have in this model. Note also that both 
bulk mass m# and localised mass m^ suffer from quadratically divergent loop corrections: 
therefore, in UED models, the naturalness problem is doublet compared to the 4D SM. 

For the (n, 0) mode, the Lagrangian contains the following mass term, including the 
contribution of the VEV: 

\v 2 \ 

i 2 m 2 KK (l - 8+) + m\ + 2m\ + — J 0{„ )O) 0( n ,o) - Aw 2 (Re^^) 2 , (3.63) 
where 0° is the down component of the doublet and 

contains the log-divergent contributions. We can now expand the Higgs doublet as 

M) = ( H(n,o) ( +i°j, 0) ) ( 3 - 65 ) 
whose components have mass 

m ^(n,0) = ^^A 1 ~ fy) + m b > m H {n ,0) = m ^(n,0) + ^ ■ ( 3 ' 66 ) 

The fields 7r ± ' are the would-be goldstone bosons that would provide the longitudinal polari- 
sation to the massive gauge bosons if the Higgs VEV were the only source of mass. However, 
for the (n, 0) modes, the main contributions to the mass comes from the extra dimension, 
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and the role of the goldstone bosons is played by the extra polarisations of the vectors. The 
Higgs VEV will in general mix the 7r =l= ' with the gauge scalars: one effect of such mixing is 
to introduce a new correction to the mass of the physical states s^ n ° Q y The correct masses 
are therefore given by: 

m 2 ± = ml, , + m%j , m\ = mi s + m\ . (3.67) 

S f n ,0) W S (n,0) *(«.<>) Z V ' 

For the (0, n) modes, the same formulas apply with the simple substitution m KK — > ^m KK . 

The localised counter-terms will also generate a mixing between states with different KK 
number. An important consequence is that the VEV spreads to other KK states and not only 
to the zero mode. Such corrections are can be relevant, as they will generate mixing terms 
between zero modes and heavy states, thus generating corrections to electroweak observables 
at tree level. Among others, mixing between the (n, 0) and (0, n) states will be generated. 
As the quartic coupling is in the bulk, the off-diagonal terms will be the same for the 4 Higgs 
states: 

5m l,(n,o)-{o,n) = ^n 2 ™M^ + 2m b ■ (3-68) 



£ dependence 

The masses of the KK Higgs bosons are only sensitive to the localised divergent contribution 
of the loops, which are proportional to the masses of the modes. A very mild £ dependence 
remains via the small finite and n-dependent terms appearing with the glides, however such 
terms are negligibly small. The only term that will introduce a £ dependence in the relative 
mass contribution is the localised mass m&, which is a constant: this contribution can be 
dealt with on the same footing as the Higgs VEV one, and not like a loop contribution. 



Tier mixing 

The structure of the tier mixing for the Higgs bosons is similar to the gauge boson ones, 
therefore he same general formulas can be used in this case. 

In the degenerate radii case £ = 1, the diagonal masses are the same, therefore as before 
the mass eigenvalues are given by sum and difference with eigenvalues 

m l,n+ = m l,(n,0) + 6m l,(n,0)-(0,n) = ^^kA 1 ~ + 3mjj , (3.69) 

m l,n- = m l,(n,0) - Sm l,(n,0)-(0,n) = ^^KK ~ m b ■ (3-70) 

As before, m Iin . is insensitive to the logarithmic divergences, however it depends on the 
localised mass mj. The contribution of the Higgs VEV to the masses of H, s° and s^, are 
given by the same formulas as in the non- degenerate case. 



4 Numerical spectra and Monte Carlo implementation 

We now turn our attention specifically to a UED incarnation of the Real Projective plane, 
where a 6D field is associated to each SM field and the SM Lagrangian is extended to 6D. 



24 



Both tiers (2, 0) and (0, 2) contain a massive vector boson for each SM gauge boson, a Dirac 
fermion for each chiral SM fermion and 4 scalars (H, S° and S ± ) corresponding to the Higgs 
field. 

Neglecting higher order operators, the model has 4 free parameters: the mass scale 
mKK = I/-R4, the asymmetry factor £ = R4/R5, the cut-off A and the localised Higgs 
mass m 2 ,. The cut-off enters logarithmically the loop corrections to masses and the loop 
induced KK momentum violating couplings. Naive dimensional analysis suggests that A ~ 
5 -i- 10 ■ rriKKi an d in the following we will always fix A = 10 m^K, a situation in which 
the mass splittings are maximal. The localised Higgs mass contribution m 2 only enters the 
mass of the Higgs resonances which are not very relevant for the LHC phenomenology being 
rarely produced and difficult to detect. For simplicity, we will fix mj = GeV from now 
on (note that m 2 may be negative and drive the electroweak symmetry breaking in the SM 
Higgs sector). 

The physical spectrum of the even tiers, (2,0) and (0,2), depends in a non trivial way 
on the ratio of the radii £. Furthermore, when the difference between the two radii is of the 
same order as the loop induced corrections to the masses, significant mixing between the 
two tiers occurs. To better illustrate this behaviour we focused on a benchmark point with 
ttikk — 500 GeV: in Table [2] we show, in the first 4 columns, the mass eigenstates of the even 
tiers for 4 values of £ = 1 (symmetric case), £ = 1.1, 1.2 and 1.5. For each particle, there 
are two values which correspond to the mass eigenstates of the mixed system (2,0)-(0,2). 
For the neutral gauge bosons, there is also a mixing induced by the Higgs VEV between 
the hypercharge B and the neutral SU(2) W 3 gauge bosons, thus the 4 mass eigenstates 
of the complete system are quoted. Similarly, for the top, the Yukawa coupling induces a 
mixing between the singlet and doublet top resonances (we here neglect the effect of the 

other Yukawa couplings), thus 4 eigenvalues are quoted. To understand the spectrum, let's 

(2) 

consider first the case of the singlet charged leptons : for increasing £, we can clearly see 
the appearance of two states whose mass is proportional to uikk and C m KK respectively. 
Moreover, for large £, the relative correction is the same for the two states, i.e. 1.0038 • itikk 
and 1.0038 • ^m KK . This shows a clear separation between the two states (2,0) and (0,2), 
and also that the correction is approximately independent on £. The largest loop corrections 
arises for the gluon and it amounts to about 10% of the mass: this means that for £ > 1.1 
the two states should already be separated. In the symmetric case, the mixing is maximal 
and we observe the separation of a set of states with very small corrections and a set with 
large corrections (the latter being highlighted in bold in the table): the two sets correspond 
to the difference and sum mass eigenstates. The former is only sensitive to finite corrections, 
therefore the mass splittings are very small, while the latter is affected by twice the log 
divergent contribution. In fact, the relative correction is approximately doubled compared 
to the asymmetric case: for the e s , for instance, the correction is 0.75%. This separation is 
not a general property of the RP 2 , but, as discussed in the previous sections, it derives from 
the fact that there are no bosonic fields with parity p r = — 1 in the model. Another important 
consequence of this accidental occurrence is that the loop induced couplings between states 
in the tiers (2,0) and (0,2) to SM particles are the same. Therefore, only the sum mass 
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Table 2: Spectra of the even tiers (2,0) and (0,2) for m KK = 500 GeV: in the first 4 columns the exact 
results including one loop and Higgs VEV contributions; in the latter 2 columns the approximate results 
implemented in CalcHEP for the degenerate (symmetric) case and asymmetric one. In bold, we highlighted 
the 2 + = (2, 0) + (0, 2) combination in the case £ = 1, and the states whose mass is proportional to ttlkk in 
the other cases: the bold numbers are the ones to be compared with the CalcHEP values. For the Higgses, 
we set m;, = an d m H — 125 GeV. 



eigenstates (the one with the largest mass splitting) couples to SM particles at one loop 
level. In other words, only the states with large splittings will be able to decay into a pair 
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of SM states. The states with small splittings can decay neither into SM states at one loop 
level, nor into a pair of odd states due to the small mass corrections. This is important for 
the phenomenology, because this kind of states will only be allowed to chain decay into the 
lightest state of their kind, emitting very soft SM particles. The lightest state, which is a 
neutral vector boson, is not protected by any symmetry thus it will eventually decay into SM 
states via higher order operators. This may lead to interesting resonant final states, however 
the branching rations and lifetime will depend on unknown parameters. For the time being, 
we can safely assume that the lifetime of the neutral vector boson is large enough for this 
particle to escape the detector without leaving any signatures. Furthermore, such states will 
not play any significant role in the phenomenology of the odd states, as they cannot decay 
into them and they cannot participate as s-channel resonance in their pair production. 
From this discussion, we can therefore identify three distinct physical regimes: 

- decoupling regime (asymmetric radii): when £ is above one, the two sets (n, 0) and 
(0,n) decouple, with the latter being heavier. The mass corrections, relative to the 
tree level mass of the tier, are approximately the same. 

- degenerate regime (symmetric radii): when £ = 1, in the UED model, the (2,0)- 
(0, 2) system splits into two subsets, of which the one with larger loop corrections 
is phenomenologically interesting. The odd states are exactly degenerate and have 
identical couplings. 

- intermediate case: in this case, the mixing is not trivial. However, it only happens 
when the radii are different but close enough, therefore in a very small region of the 
parameter space. 



4.1 Monte Carlo implementation 

We implemented this model in the FeynRules [16] package of Mathematica, which allows 
to interface the model with Monte Carlo generators like CalcHEP, MadGraph4 and 5. The 
implementation should contain the full tiers (1,0), (0,1), (2,0), (0,2) with ^-dependent 
mass corrections and mixing between the two even tiers. The tier (1, 1) may also give rise 
to interesting phenomenology, however there is a dependence on the unknown coefficients of 
the localised Lagrangians, which allow for the direct decays into a SM pair: an example of a 



spectacular 4-top signal has been studied in 17 . In the following we will not consider this 



tier any further. The full implementation would contain a huge number of new states and 
interactions. Following the above discussion of the three physical regimes, however, we can 
simplify the model by considering two separate limits: the case of asymmetric radii and the 
degenerate case. 



Asymmetric regime 

For asymmetric radii, £ > 1.1, the two sets (n, 0) and (0, n) decouple because the mixing 
between the even states are negligible. Furthermore, the relative mass corrections are nearly 
independent on £: the model reduces to two uncorrelated sets, characterised my tukk = 1/R 
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and rriKK — £/R- To simplify the implementation, therefore, we can consider a model with 
only one set of odd (1, 0) states and one set of even (2, 0) modes (for convenience, we can 
label them with a (1) and (2) dropping the zero). The heavier states (0, 1) and (0, 2) can be 
described by the same simplified model with rescaled mass scale m KK . 

Thus, in order to study the phenomenology of a model with uikk and £, it suffices to 
generate events with mass scales uikk and £ uikk and sum the two results. 

Degenerate (symmetric) regime 

In the symmetric case, the model contains two degenerate and identical odd tiers, while 
the even tiers consist of one (the sum) that couples to SM states and has larger mass split- 
tings, and one (the difference) that does not. The latter is not relevant for the phenomenology 
we are interested in. To simplify the implementation, we can study a model that includes one 
of the two odd states, and the even states with larger mass splittings: the particle content is 
therefore the same as in the asymmetric simplified model! In the mass corrections, we need 
to supplement a factor of 2 in front of the log divergent contributions to the even masses 
only. Furthermore, when studying the phenomenology of the odd tier, we need to multiply 
the results by a factor of 2 to take into account the presence of two degenerate and identical 
tiers. The couplings of the odd tier will also be affected by the change of basis: here we 
are interested uniquely in the phenomenology of even and odd tiers. The processes we will 
consider are: 

- pair production or annihilation of odd tiers, which are sensitive to the even tier via the 
couplings (2, 0)(1, 0)(1, 0) and (0, 2)(0, 1)(0, 1) in diagrams with a resonant s-channel 
even state. In the mass eigenstate, the two couplings pick up a factor of l/y/2. How- 
ever, this coupling will also enter the decays of the even state into a pair of odd states: 
in our simplified model which includes only one odd state, this partial width would 
be underestimated by a factor of two. To compensate for this, and have the correct 
branching ratios, we do not include the 1 / y/2 in the coupling. Therefore, when com- 
puting resonant production of the odd states, or annihilation via a resonant even state, 
it is intended that the sum over the two degenerate odd tiers is implicitly done. 

- pair production of the even states are sensitive to the couplings (0, 0)(2, 0)(2, 0) and 
(0, 0)(0, 2)(0, 2), which are equal at three level. In the mass eigenstates, no extra factor 
is generated. These couplings will also generate decays (2) — > (2)' + SM. 

- single production of the even state: this process is sensitive to the loop induced cou- 
plings (2,0) (0,0) (0,0) and (0, 2)(0, 0)(0, 0), which have the same coefficient. In the 
mass eigenstate, the coupling has an extra factor y/2. This coupling will also enter the 
decays (2) ->■ SM. 

We have implemented the two simplified models in FeynRules, and interfaced them with 
both MadGraph and CalcHEP. As the field content and couplings are the same, we have 
implemented a single model and introduced a parameter that distinguishes the two cases 
by adding the appropriate a/2 factors as described above. In the implementation of the 
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mass corrections, we simply use the formulas for £ = 1, which are formally valid only in the 
symmetric case, relying on the fact the mass corrections show a very mild dependence on 
£. In the last two columns in Table |2j we show the spectra as computed in the CalcHEP 
implementation of the model. We see that the approximate masses closely reproduce both 
the asymmetric and symmetric cases. The model also contains the full set of loop induced 
couplings, as computed in Ref. |5], and the full one loop masses for both even and odd states. 

5 Estimates of LHC 2011 bounds 

The LHC experiments have collected an impressive amount of data, 5 fb _1 per experiment [j 
in 2011 at a centre of mass energy of 7 TeV and a similar amount so far at 8 TeV in 2012 
with a forecast integrated luminosity of 30 fb _1 collected by the end of 2012, before the 
2 year technical stop. Even though no search specific to this model have been performed, 
many of the searches designed for supersymmetry or other exotic signals and models can be 
used to pose severe bounds on the RP 2 . 

The lightest odd tiers will give rise to signatures similar to supersymmetry due to the 
decay chains ending inevitably into a stable neutral particle (missing transverse momentum). 
However, the spectrum of the tiers is rather compressed due to the smallness of the splittings: 
the most effective searches should be the ones looking for supersymmetric signals in events 
with jets and missing transverse energy where, in this case, the jets are coming mainly from 
initial state radiation. A detailed simulation is essential in order to estimate the efficiency 
of the experimental cuts on the signal events, thus a study specific to this model is in 
progress [fj. 

In the following we will focus on the phenomenology of the even tiers (2,0) and (0,2): 
the main interest here is that, due to loop induced vertices, such states can decay directly 
into a pair of SM particles, thus giving rise to resonant signatures. Even though the mass 
of the even tiers is roughly double compared to the lightest odd tiers, thus the productions 
cross sections are considerably smaller, the strong bounds on resonances can be comparable 
or even stronger than the bounds from missing energy searches. In this model, only bosons 
can couple directly to a pair of SM states, giving rise to potential resonances in two fermions 
or two gauge bosons or Higgses. Fermions can potentially couple to a SM Higgs (or massive 
gauge boson) and a SM fermion, however such coupling will be proportional to the Yukawa 
coupling, thus it is numerically relevant only for third generation quarks. Furthermore, 
resonant final states with a Higgs or massive gauge boson plus a quark are very difficult to 
reconstruct, therefore we will ignore this possibility in the following. The interesting final 
states are: 

8 Here we focus on the general purpose ATLAS and CMS detectors. 
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There exists searches looking for resonances in all these final states, however the most sensi- 
tive ones are the leptonic searches for Z' — > and W — > l + v resonances. This holds true 
in this model too, even though the branching ratios into leptons are very small. Jets, which 
have much larger recurrence, need to fight against a very large QCD background: searches of 
dijet resonances are only sensitive to rather large invariant masses (above 1 TeV) and order 
pb cross sections. On the other hand, there is an interesting CMS search for a pair of dijet 
resonances which will give interesting results in this case. Resonances involving gauge bosons 
or tops pay the price of leptonic branching ratios for the W and/or Z decays, and severe cuts 
to reconstruct the heavy states, thus they will not be competitive. In the following, we will 
try to estimate the reach of the relevant searches without performing a complete simulation 
but just counting the effective cross section for final states with the resonance and assuming 
the same experimental acceptance as for the model analysed in the experimental paper. The 
limitations of this approach will be discussed case by case. 

5.1 Simulation details 

For the numerical study we make use of the Madgraph and CalcHEP interface of the model 
via FeynRules, described in the previous section. All the cross sections used here were 
computed using Madgraph, and checked with CalcHEP. We set the two free parameters of 
the model to AR = 10 and m b = GeV: the former controls the loop induced mass splitting 
and we fix it to the maximum value we would expect from naive dimensional analysis in 
6D, the latter controls the mass of the Higgs resonances. The Higgs resonances turn out 
to be irrelevant for the channels we are interested in here, thus changing mj will not affect 
our results significantly. In all the calculations we consider both the symmetric radii case 
(i?4 = -R5, or £ = 1) and the asymmetric radii case (i? 4 ^> R 5 , or £ ^> 1), as in the simplified 
FeynRules implementation. 

The resonant channels only come from electroweak gauge bosons -4(2), Z@) and or the 
gluon G{2)- The latter is the heaviest state in the tier, thus it can only be directly produced, 
while the electroweak bosons can derive form the decays of heavier states (the quarks). 
Direct pair production of electroweak bosons is largely overwhelmed by the production via 
decays of coloured states, therefore we will only focus on the latter and neglect subdominant 
electroweak processes. On the other hand, it is possible to singly produce all the gauge 
bosons via the loop induced coupling to the quarks: the resonant states can thus come from 
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Figure 5: Pair production cross sections of even heavy states for asymmetric (left two plots) and symmetric 
(right two plots) radii. The colours correspond to: total (black), QQ (solid green), GQ (magenta), GG (red), 
and QQ3 (dashed green). The solid green line is barely visible, as it almost coincides with the black one. 



resonant Drell Yan production or from a chain decay of the gauge boson. We will consider 
the following production channels: 

- Pair production: 

a) pair of first and second generation quarks QQ , including the processes (?(2)(?(2)> 
<7(2)<?(2) and q(2)Q(2) with all the flavour combinations; 

b) production of a quark in association with a gluon GQ, where the heavy quark is 
a first or second generation partner; 

c) pair production of gluons GG; 

d) pair production of third generation quarks QQ3, including &(2)&(2) and £(2)?(2)- 

- Single production: 

a) single production of a gluon G(2), which can decay into leptons via a chain decay 
involving heavy quarks; 

b) single production of Z^ 2 ) and W^, which can decay into leptons (or jets) either 
directly or via heavy leptons; 

c) single production of a ^4(2) > which can only decay directly into a pair of leptons 
(or jets). 
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Figure 6: Single production cross sections of even heavy states for asymmetric (left two plots) and sym- 
metric (right two plots) radii. The colours correspond to: gluon Gn) (red), Wn) (green), Zr$\ (blue) and 
A( 2 ) (grey). For comparison, the dotted line is the total pair cross section. 

All those channels can produce resonant SM particle pairs at the end of a short or long decay 
chain. Some examples of decay chains, with a resonant / u + /i~ pair, are: 

* gg } ->• q'W^ + ->• gVeg )+ q'ue+A® g Ve+ [pL + pL~); 

* ->• qq ( n ] -> ggZ( 2 ) ->• gg[/iV~]; 

* _). /i+z/g^ -> /i+/i-A( 2 ) -> [//+//-]; 

where the two resonant leptons are in brackets. The decay chain can contain quite a few 
additional SM particles, including jets and leptons, which can in principle complicate the 
reconstruction of the resonant pair. However, the small mass difference between states in 
the same tier will ensure that the momentum carried by the intermediate products in the 
chain is much smaller than the momenta of the resonant pair, which inherits the mass of 
the heavy gauge boson. Therefore, it is a reasonable assumption that the soft intermediate 
decay products will not affect significantly the acceptance of the experimental cuts. We will 
base our analysis on this assumption, leaving a detailed simulation for future work. We also 
limit ourselves to leading order cross sections, even though QCD corrections may give rise to 
significant increase (or decrease) of the cross sections. The results for the total cross sections 



32 



500 

ma 



500 

»>KK 



Figure 7: Ratio of the 8 TeV to the 7 TeV cross section for pair production (left) and single production 
(right) channels. The same color codes as in the previous figures is adopted. The shown plots refer to the 
asymmetric case, similar values are obtained in the symmetric one. 

are shown in Figure [5] and [6j From Figure [5] we can see that all pair production cross sections 
are dominated by the QQ channel, and that they drop very quickly with increasing tukk' 
this is due to the fact that two states of mass ~ 2 * m KK must be produced. Furthermore, 
the symmetric radii case fares slightly smaller cross sections, due to the fact that the masses 
are larger than the asymmetric case because of larger loop contributions. We also see a 
significant increase from 7 TeV centre of mass energy to 8 TeV. From Figure |6j on the 
other hand, we can see that single production cross sections decrease much more slowly with 
increasing m^K, until they dominate for large masses {vtikk ^ 450-^500 GeV for asymmetric 
radii). In the symmetric radii case, the single production cross section is larger: this is due 
to the larger loop induced coupling. Single production, therefore, tends to dominate over 
pair production for smaller radii {itlkk ^ 350 -j- 400 GeV). We also see that the increase of 
cross section to 8 TeV is not as significant as in the pair production case. This fact is more 
clear in Figure [7| where we plot the gain from 7 to 8 TeV in the asymmetric case. For pair 
production, the gain can be significantly larger than 1 and is particularly effective for final 
states involving gluons, while for single production channels the gain is always limited. 

The main difference between the asymmetric and symmetric radii cases so far is that in 
the symmetric case the single production cross sections are much more important. Another 
fundamental difference appears in the decay rates: in the asymmetric case, decays into odd 
gluons are kinematically forbidden because the mass correction to the odd gluon is too large. 
In the symmetric case, the mass corrections to even (2) states are doubled. As a consequence, 
decays into odd gluons open up: this means that for even quarks and gluons, the channel 

— > G^q^ and — > G^G^ will significantly increase the decay rate into odd states. 
Therefore, rates into resonant pairs will be significantly reduced, particularly so in channels 
involving heavy quarks. 
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Figure 8: Leptonic Z' bounds in the asymmetric (left) and symmetric (right) cases. The solid black line 
represents the total cross section, while the dashed (dash-dotted) one shows the pair (single) production 
component. The blue line is, for comparison, the 8 TeV total cross sections. The CMS (ATLAS) excluded 
region is shown in yellow (green). 

5.2 Z' — > l + l~ searches 

The most promising search is for a resonance in the dilepton channel, which is one of the 



cleanest channels at hadronic colliders like the LHC. Both CMS 18 and ATLAS 19 have 



published the final analysis of the full 2011 dataset at 7 TeV, while only CMS has published 



preliminary results with the 8 TeV run data 20 . Both experiments search for final states 



with at least two energetic electrons or muons with opposite sign, and reconstruct the in- 
variant mass of the pair. The electron channel is slightly more efficient than the muon one, 
because the total energy can be measured in the electromagnetic calorimeter. On the other 
hand, the muon momentum can only be reconstructed by measuring the curvature of the 
track in the muon chambers, and this procedure becomes more difficult for energetic muons. 
The expected signal is a Gaussian bump in the dilepton invariant mass distribution, and the 
efficiency of the signal is estimated by simulating a narrow resonance for which the width of 
the Gaussian peak is determined by the experimental resolution on the lepton momentum 
reconstruction. 

In the model under investigation, the resonant dilepton pair is accompanied by extra soft 
leptons coming from the intermediate decay chain, however, the momenta of the resonant 
pair should be much larger than the soft intermediate ones. Therefore, the presence of extra 
leptons should not affect the reconstruction of the peak. Furthermore, the decaying bosons 
are very narrow, with widths of the order of 0.005% of the mass for and 0.01% for 
the Z( 2 \ The most serious issue is that the resonant lepton pairs come from two different 
states with different masses: the loop corrections generate a mass splitting between the 
and the which amounts to about 5% of ttikk (~ 2.5% of the mass) in the asymmetric 
case and 9.5% of tukk (~ 5% of the mass) in the symmetric case. The CMS search at 7 
TeV 1 18 1 reports that the resolution on the dimuon invariant mass is 6.5% at masses around 
1 TeV, increasing to 12% at 2 TeV, while the resolution in the dielectron channel is 1.1% 
for electron in the barrel and 2.3% when one electron is in the endcap calorimeter, staying 
constant above 500 GeV. These numbers show that the electron channel has enough resolu- 
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tion to disentangle the two resonances, and a dedicated search may give very tight bounds 
on this model. For the muon channel, the experimental resolution is still too large. The 
experimental bounds on a Z' mass is performed by using an unbinned likelihood analysis, 
where the signal is modelled by a Breit-Wigner function convoluted with a Gaussian resolu- 
tion function. The width in the Breit-Wigner function is set to be equal to a sequential Z', 
thus amounting to 3.1% of the mass. The single bump used in setting the bounds would in 
all cases encompass both resonances in this model, therefore in the following we will assume 
that both resonances contribute to the same bump and sum the two cross sections when 
comparing the expected cross section to the experimental bound. A similar strategy is used 



by the ATLAS collaboration 19 : for their detector, the energy resolution on electrons is 
estimated to be 1.2% in the barrel and 1.8% in the endcap regions, while for muons the pt 
resolution ranges from 10% to 25%. 

In this model, the effective branching ratios into a resonant pair of leptons are very small: 
they range from 0.2% to 0.5% for coloured states in the asymmetric radii case, and from 
0.05% to 0.3% in the symmetric case. Nevertheless, the total cross section is still quite large 
and comparable to the experimental bounds. The results of our calculation are shown in 
Figure |8j where the black lines show the total cross section while the dashed (dot-dashed) 
ones show the component from pair (single) production channels only. The experimental 
bounds, in green for ATLAS and yellow for CMS, refer to the combined results with 5 fb _1 
of data at 7 TeV, where we can see that at large masses the CMS bound prevails. In the 
asymmetric case, the best bound is at vhkk ~ 575 GeV, which is in a mass region where the 
pair and single production channels are comparable. In the symmetric case (right plot), we 
see that the pair production channel is highly suppressed because of the heavier mass, while 
the single channel is suppressed to a smaller degree due to the smaller branching ratios. The 
net result is a bound at around m KK ~ 440 GeV, dominated by pair production. 

The CMS collaboration has also published preliminary results on an analysis performed 
on the 2012 data at 8 TeV, however no combined electron-muon results is reported. The 
best quoted bound comes from a combination of the 7 and 8 TeV datasets, which is done 
by rescaling the signal at 7 TeV to 8 TeV. As this rescaling is model dependent, we cannot 
directly use this result here. Nevertheless, in the Figure [8j we show for comparison the 
expected cross section at 8 TeV in blue 

5.3 W — > Iv searches 

The second most promising analysis is the search for a W in the leptonic channel. Both 



CMS 21 and ATLAS 22 collaborations have published their final searches with the 2011 



7 TeV dataset, while only CMS 23 has preliminary results available for the 2012 8 TeV 
dataset. This search is based on the requirement of a single high transverse momentum 
electron or muon in association with missing transverse momentum (which may derive from 
the neutrino). The main observable used to isolate the signal is the transverse mass Mt, 
defined as 



M T = \ 2 ■ p l T ■ p^ lss ■ (1 - cos A<j)) , (5.71) 
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Figure 9: Leptonic W bounds in the asymmetric (left) and symmetric (right) cases. The solid lines 
represent the total cross section, while the dashed (dash-dotted) one shows the pair (single) production 
component. The dotted line is the total cross section without the MET veto. The upper plots refer to the 
7 TeV data, while the lower ones to the 8 TeV CMS analysis. The CMS (ATLAS) excluded region is shown 
in yellow (green). 

where p l T is the transverse momentum of the electron, p™ lss is the missing transverse momen- 
tum in the event and S(p is the azimuthal angle between the lepton and missing transverse 
momenta. A signal from a massive W — > Iv would produce a Jacobian peak in the Mt 
distribution at the W mass, therefore a lower cut on Mt is imposed somewhat below the 
W mass used in the simulated signal. As before, the presence of intermediate decay prod- 
ucts should not significantly affect this search. However, in the pair production channels, 
the decays of the second heavy state may contain additional invisible particles which will 
contribute to the total p™ lss of the event and distort the Mt distribution. The most worri- 
some cases involve the decay to a pair of odd states, which have significantly large branching 
rations. To be conservative, in this analysis we will veto on the decay of the second heavy 
state into an odd pair, nevertheless the result without this veto are also shown. The actual 
bound should lay in between these two lines. 

Both CMS and ATLAS extract the bounds on a W — > Iv that does not interfere with 
the SM, and for each mass point analysed they impose an optimised cut on Mt- Here for 
simplicity we will assume that the efficiency of the experimental cuts are the same as the ones 
in the model they consider, and we simply compare the value of our effective cross section 
to the bound. The results are shown in Figure [9] for both asymmetric and symmetric radii, 
and for 7 and 8 TeV studies. At 7 TeV, the bounds from ATLAS and CMS are comparable. 
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In the asymmetric case, the bound is significantly lower than the the dilepton channels, and 
it amounts to m^K ~ 430 (435) GeV for the 7 (8) TeV analysis. The bounds is in a region 
dominated by the pair production, and the inclusion of events with a pair of odd states 
would improve the bound by about 50 GeV. In the asymmetric before, we can see 

a significant decrease of the pair production channels and a less important decrease in the 
single production. The decrease is however sufficient to remove the bound, as the predicted 
cross section always lies below the excluded region. The inclusion of events with an odd pair 
would significantly enhance the event rate, with a maximum bound of tukk ~ 300 GeV. 

5.4 Paired dijet search 

The coloured heavy states have rather large branching ratios in events with a resonant jet 
pair, ranging from 35% to 65% in the asymmetric case, and from 12% to 23% in the symmetric 
case. However, analyses based on jets suffer from large QCD background which are difficult 
to deal with. Searches for dijet resonances have been published from both ATLAS and CMS, 
however the bound on the cross section is significant only in regions with very high invariant 
mass (above a TeV). On the other hand, CMS has published an interesting search of events 
with a pair of dijet resonances with similar mass [24], based on 2.2 fb _1 integrated luminosity 
from the 2011 data. Due to the large branching ratios, the pair production channels are due 
to produce a fairly large number of such states, therefore this search may be competitive for 
this model. 

The search strategy is the following: the selected events must have at least 4 jets with 
transverse momentum p? > 150 GeV. The four highest px jets are combined into two pairs 
with a large separation between the two pairs in order to minimise overlaps, and, out of all 
the combinations, the search selects the one with the largest Am/m avg , where Am is the 
difference in invariant mass and m avg is the average invariant mass of the two pairs. Finally, 
only events with Am/m avg < 15% are selected in order to minimise the QCD multijet 
background. An additional cut on 5 = PT, pa iri + PT, P air2 — mavg > 25 GeV is also imposed 
to further reduce the QCD background and have a smooth tail of background events at 
n^avg > 320 GeV. These selections do not suffer from the presence of the soft intermediate 
products, however a significant factor in the game is the fact that the dijet pairs in this model 
may arise from 4 states with different masses: A^ 2 \ Z^ 2 \ and G^ 2 \ The mass slitting 
between the and with the lighter is of the order of 2.5% (5%) of the mass in the 
asymmetric (symmetric) case, while the splitting of the gluon with is of the order 
of 10% (18%). While the splitting between electroweak gauge bosons is small compared to 
the 15% cut on Am/m avg , events with gluons may fall outside of the selected region. This 
problem is not severe, though, because the signal is dominated by the electroweak gauge 
bosons: more than 80% of the signal, therefore the signal loss would be limited. The reason 
for this is that the G^ is the heaviest particle in the tier, therefore it cannot appear in 
decay chains of lighter states like quarks. Thus, the signal with a G^ resonance paired with 
an electroweak boson can only appear in the GQ or GG channels, with one gluon decaying 
directly into a pair of jets. In the following we will neglect this potential signal loss, and 
consider the total cross section. 
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Figure 10: Di-jet pair bounds in the asymmetric (left) and symmetric (right) cases. The solid black lines 
show the total effective cross section assuming an experimental efficiency of 3% and 13% respectively. The 
dashed blue lines are, for comparison, the 8 TeV total cross sections. The CMS excluded region is shown in 
yellow. 



Another important issue in this channel is the fact that the published bound is on the 
production cross section times branching ratio times the experimental acceptance. The 
model considered in the CMS search is a model with scalar colour octets (coloron), which 
can decay into a pair of quarks qq, and they are pair produced via QCD interactions. The 
experimental acceptance for this model is 3% for a coloron mass of 300 GeV, raising to 13% 
for a mass of 1 TeV. The search has been performed for masses in the range 320 to 1200 
GeV. To extract a reliable bound from this search, a detailed simulation of the signal is 
indispensable, in order to extract the experimental acceptance. However, in order to have 
an estimate of the bound, we will here assume that the acceptance is within 3% to 13%, and 
leave a detailed simulation for future work. 

The result of our calculations are shown in Figure 10 where the two black lines correspond 
to the two coloron acceptances, and the CMS excluded region in yellow is based on a 7 TeV 
dataset of 2.2 fb _1 . The curves show that in the asymmetric case, this search is potentially 
relevant due to the fairly large cross sections obtained in this case. Assuming the coloron 
acceptance, the bound on txikk is expected to range from 380 to 540 GeV, a region which is 
close to the dilepton channel. Moreover, this bound is only based on half of the total 2011 
dataset. The blue dashed lines in the plot show the predicted cross section at 8 TeV. In 
the symmetric case, the prospect is much more difficult: in fact, the cross section is heavily 
suppressed because it only receives contribution from the pair production channels, and we 
expect a bound in the 300 GeV mass region, thus well below the dilepton channel. 



6 Conclusions 

The excellent performance of the LHC and of the related experiments has allowed to collect a 
remarkable amount of high quality data, 5 fb _1 at 7 TeV in 2011 and an expected 30 fb _1 at 8 
TeV by the end of 2012. This has allowed for the discovery of a new state that resembles the 
Standard Model Higgs boson, and to place very strong bounds on the masses of constrained 
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Channel 


asymmetric 


symmetric 


Z' ->■ i+r 

W ->■ Zz/ 
paired dijet 


575 
430 
380 -=- 540 


440 
270 -r 295 



Table 3: Estimated bounds on ttikk in the asymmetric (£ 3> 1) and symmetric (£ = 1) cases. 

supersymmetric models and other simple extensions of the SM. The experimental results, 
which does not show any sign of physics beyond the predictions of the SM, however, are yet 
not sufficient to disfavour the naturalness argument which would prefer new states at the 
TeV scale. In fact, more work and data are necessary to strongly constraint top partners or 
alternative models to supersymmetry. In this paper we focused on one alternative scenario, 
where the Dark Matter particle that fills the Universe is provided by extra dimensions. In 
fact, symmetries of the compact space, relic of the extended Poincare invariance after the 
compactification, can ensure the stability of the lightest Kaluza-Klein resonance of the SM 
states. The presence of such a symmetry is indeed a remarkably strong requirement, allowing 
for instance to select a single flat orbifold in 5 or 6 dimensional models: the Real Projective 
plane (RP 2 ). 

In this paper we study a UED realisation of the RP 2 , where a 6D field is associated to each 
SM field, including the Higgs field. Loop corrections are essential to study the phenomenology 
of this class of models. In this work, we completed the one loop analysis of the first odd and 
even tiers, by presenting the complete calculation of the one-loop corrections to the masses 
of the even tiers, and we generalised the results for the odd tiers [2] to the case of asymmetric 
radii. The results for the loop induced couplings of the even tiers can be found in [5] . These 
loop results are essential to study the phenomenology of the model, as they determine the 
decay and production rates of the new states. In this work, we estimated the bounds from 
resonant decays of even gauge bosons into a pair of standard model particles. We identified 
3 promising channels: dilepton resonances, decays into a single lepton and neutrino, and 
decay into dijet of a pair of states. In the estimate, we simply count the effective cross 
section of the final states of interest, and we assume that the experimental efficiency is the 
same as the model used by the experimental groups to publish the bounds. This procedure is 
approximate, and a detailed simulation would be necessary to extract more reliable numbers; 
nevertheless, the estimates can give us an idea of the relevance of the bounds. The results are 
summarised in the Table |3j the most promising channel is the dilepton resonance, that gives 
a lower bound on m KK of 575 (440) GeV in the asymmetric (degenerate) case. However, 
to extract the estimate we assumed that all the events fall under the same Gaussian peak, 
while in reality they arise from two separate peaks from Zi 2 ) and Ann- A fit to the data 
that takes into account the presence of two peaks may improve the bound, especially in the 
electron channel that has better energy resolution. A surprisingly competitive channel is the 
search for paired dijet resonances: in this case, the bound depends crucially on the efficiency 
of the selection rules. Furthermore, the CMS study is only based on 2.2 fb _1 of data at 7 
TeV, thus the bound may be significantly improved if the QCD background is properly taken 
into account. Our estimate, based on the efficiencies quoted on the CMS paper and based 
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on a scalar coloron model, show that in both cases this search may be competitive with the 
dilepton one. Other channels, like ti resonances, tb resonances, WV ones and dijet, are not 
competitive as they are mostly sensitive to large invariant masses and large cross section. 
The price to be paid here is either a small leptonic branching ratio of the W and Z, or the 
large QCD background. The resonant channels are very important because they simplify the 
searches, and also because such signatures are not present in supersymmetric models. The 
UED model also have more supersymmetry-like signatures coming from the chain decays of 
the odd tiers, ending into a stable neutral particle. In this case, the signal will resemble that 
of a compressed supersymmetric spectrum, and the most sensitive channel should be the one 
with jets and missing transverse momentum. A detailed analysis of the reach of this channel 
is in progress 

Contrary to supersymmetry, the masses of the KK resonances cannot be pushed to an 
arbitrary large mass scale: in fact, the calculation of the relic abundance selects a finite 
range for the preferred value of m^K, and the overclosure of the Universe sets an upper 
bound. Even considering the uncertainties related to our knowledge of the Cosmological 
model, this feature persists. In Ref. [2j the preferred range has been estimated using an 
analytical approximation and only considering co-annihilation of all the states in the odd 
tier: the result is a preferred range of 300 < m KK < 400 (210 < m KK < 280) for asymmetric 
(degenerate) radii. However, this rough estimate is not enough, as it does neglect the impact 
of the even tiers 25 , that tend to push the preferred range towards higher values. A detailed 



analysis of these effects will be presented in a forthcoming publication 10 . 

The RP 2 also offer other phenomenologically interesting aspects, that have not been 
explored yet. For instance, the odd tiers (2,1) and (1,2) have masses close to the even 
ones under study, fn^ 2 ,i) = V^i^kk ~ 2.25 m KK , thus their production rates will be similar. 
Furthermore, they contain twice as many fermions. The only allowed decay modes are decays 
within the tier, or decays to a lighter odd state plus a SM one, (2, 1) — > (0, 1) + (0,0), via 
loop induced couplings. These tiers will therefore be a novel source of events with jets and 
missing transverse momentum, and they may significantly affect the bound on the radius. 
There are also tiers that can only decay via localised counterterms, like (1, 1) or the (2)_ 
in the degenerate case: they may give rise to resonances or new sources of missing energy 
depending on the localised interactions. 
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A Appendix: loop functions 

The contribution of the torus to the loops is proportional to the function 

* m,T»2 c z y 1 ■» y 

{«1.»2) # {0,0} 

' 2 ..f^+cothf^-lV— (A.2) 
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The approximation is very good for moderately large values of £, i.e. for £ < 5. Note that 

en(l/0 = T 6 (0 , (A.3) 

thus the torus contribution to (n, 0) and (0, n) modes is the same. T 6 (£) is a rapidly increasing 
function of £: for instance, T 6 (l) ~ 1.92 while T 6 (2) ~ 6.7. 

The KK number dependent functions generated by the glides, can be expressed in terms 
of the following integrals: 

f°° k 3 

$i(n) = 2tt 3 / ^ , (A.4) 

Jo yk 2 + n 2 smh.Tiyk 2 + n 2 

$ 2 (n) = 2tt 3 /""rfjfc MV^^) (A5) 
Jo Vk 2 + n 2 sinh ir\/k 2 + n 2 

$ 3 ( n ) = 2vr 3 / dA; — v y . (A.6) 

Jo n\/k 2 + n 2 smh 7C\/k 2 + n 2 

These integrals are exponentially suppressed for large n, thus the impact is reduced for higher 
modes. 

The functions appearing in loops for gauge bosons and fermions read 

S g (x) = 2$ 1 (x) - 16$ 2 (x) -8® 3 (x), (A.7) 
S.(x) = $i(x) -4$ 2 (x) -4$ 3 {x), (A.8) 

V 9 (x) = h x (x)-^ 2 (x) + ^ 3 (x), (A.9) 

V s (x) = U 1 (x) + U 2 (x) + U 3 (x), (A.10) 
F g (x) = -$!(x)-4$ 2 (x), (A.ll) 
F s (x) = ~$i(x). (A.12) 
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